Three-dimensional strong localization of matter waves by scattering from atoms in a 

lattice with a confinement-induced resonance 
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The possibility of using ultracold atoms to observe strong localization of matter waves is now the 
subject of a great interest, as undesirable decoherence and interactions can be made negligible in 
these systems. It was proposed that a static disordered potential can be realized by trapping atoms 
of a given species in randomly chosen sites of a deep 3D optical lattice with no multiple occupation. 
We analyze in detail the prospects of this scheme for observing localized states in 3D for a matter 
wave of a different atomic species that interacts with the trapped particles and that is sufficiently 
far detuned from the optical lattice to be insensitive to it. We demonstrate that at low energy a 
large number of 3D strongly localized states can be produced for the matter wave, if the effective 
scattering length describing the interaction of the matter wave with a trapped atom is of the order of 
the mean distance between the trapped particles. Such high values of the effective scattering length 
can be obtained by using a Feshbach resonance to adjust the free space inter-species scattering 
length and by taking advantage of confinement-induced resonances induced by the trapping of the 
scatterers in the lattice. 

PACS numbers: 



I. INTRODUCTION 

The recent advances in the manipulation of ultracold 
gases have made it possible to employ these systems to 
accurately simulate several non-trivial problems in con- 
densed matter physics. As examples, we may mention the 
exploration of the BCS to BEC crossover [J and the su- 
perfluid to Mott-insulator transition Q. Disorder plays 
an important role in the theory of solid state, affecting in 
a substantial way the transport properties of various sys- 
tems. Special attention has in the past been dedicated to 
studies of light propagation in strongly-scattering pow- 
ders and electron transport in the presence of impurities 

It looks therefore interesting to introduce a controlled 
disorder in the experiments with ultracold atoms, in or- 
der to provide a closer modeling of realistic systems of 
condensed matter physics. It was indeed predicted that 
atomic gases stored in optical lattices would be good can- 
didates to experimentally observe the effect of a disor- 
dered or quasi-periodic potential on an interacting Bosc 
gas or on interacting Fermi-Bose mixture |y, Q. First 
experimental results along this research line have been 
recently reported, i.e. the observation of a Bose glass in 
a quasi-periodic potential Q and the study of spatial co- 
herence properties of an interacting Bosc gas trapped in 
a lattice in presence of a disordered ensemble of fermionic 
atoms |3- 

In this paper, we consider a variant of this line of re- 
search, that is the possibility of looking for genuine local- 
ized states of a non-interacting matter wave exposed to 
static disorder in continuous space. Localized states are 



stationary states with a square integrable wave function 
at an energy where the classical motion is not bounded 
spatially. In a paper that dates back to the early years 
of quantum mechanics, von Neumann and Wigner [lOl | 
showed that the Schrodingcr equation can admit square 
integrable eigenstates embedded in the continuum of 
states with energy higher than the maximum of the po- 
tential. After the work of Anderson [ll|, it is expected 
that disordered potentials can generically lead in 3D to a 
quantum phase transition, a macroscopic number of lo- 
calized states being present at low energy. Such a phase 
transition in 3D is not straightforward to observe, as it 
is sensitive to decoherence and wave absorption effects, 
and requires a mean free path of the wave / smaller than 
its wave length A, as stated by the loffe-Regel criterion 



A>L 



(1) 



The study of localization of light is a well developed ex- 
perimental subject: strong localization of light has been 
reported in semi-conductor powders [I^l, and weak lo- 
calization effects of light in a gas of cold atoms are the 
subject of an intensive experimental study |l3 |. On the 
contrary, for matter waves, no direct evidence of local- 
ization was obtained in 3D. Matter waves made of ul- 
tracold atoms are good candidates in this respect, due 
to their weak coupling to the environment and to the 
possibility of tuning their interactions with a Feshbach 
resonance |l4| . An open problem is however to know if 
strong enough disorder can be introduced in these gases 
to lead to reasonably short localization lengths in 3D. 

A natural way to produce a disordered potential in 
atomic gases is to use the speckle pattern of a laser beam 
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FIG. 1: (Color online) Sketch of the model of disorder con- 
sidered here: a matter wave (A, blue) scatters on randomly 
distributed B particles (red), each occupying the vibrational 
ground state of a node in a 3D optical lattice (here the av- 
erage occupancy p is 0.1). The lattice does not act on the 
matter wave A. 



numerically in Sec. IIIII that it leads to the appearance of 
a large number of localized states for a range of parame- 
ters accessible in present experiments, provided that the 
effective coupling between the matter wave and a sin- 
gle scatterer is tuned to a value of the order of the mean 
scatterer separation. Scction lTVl is dedicated to the quan- 
titative description of the scattering between the matter 
wave and a single trapped scatterer, and to the discussion 
of the confinement-induced resonances thereby arising: 
we show indeed that large enough effective coupling con- 
stants can be obtained. Our conclusions are presented 
in Sec.0 including a discussion of a possible strategy of 
observation of strong localization. 



II. OUR MODEL 



A. The disordered potential 



|l5|. Many experiments on Bose-Einstein condensates 
in ID random optical potentials have very recently been 
reported [l^ [13, lH, 113 , and they provide evidence for 
disorder-related effects such as fragmentation of the con- 
densate, suppression of diffusion, frequency shifts and 
damping of collective oscillations. Theoretically, these 
effects were discussed in |23|. Genuine strong localiza- 
tion in ID, in the non-interacting regime, has not been 
reported yet in these experiments, and the implementa- 
tion of the disordered optical potential in 3D remains to 
be done. Also the theoretical analysis of matter wave 
localization in a speckle pattern, recently performed in 
2D [21|, has not been done in a detailed way in the 3D 
strong localization regime. 

An alternative method to realize a disordered potential 
was proposed in |22|: a matter wave, made of atoms of a 
species A, scatters off a set of atoms of another species B 
that are trapped in randomly chosen sites of a deep op- 
tical lattice with no multiply occupied sites (see Fig.^l. 
As we will choose the lattice to be very far detuned for 
the species A, the matter wave moves unaffected through 
the optical lattice, and interacts only with the B atoms: 
this excludes classical localization effects in local poten- 
tial minima. The disorder can be made very strong, since 
(i) the correlation length of the disorder can be as small 
as 0.5/im (the spatial period of the optical lattice), and 
(ii) the scattering cross section of the matter wave off a 
single B atom can be made as high as allowed by quan- 
tum mechanics (the so-called unitary limit) by use of a 
A ^ B interspecies Feshbach resonance, making it pos- 
sible to dramatically reduce the mean free path of the 
matter wave. Furthermore, as we shall take advantage 
of, this model allows a straightforward exact numerical 
analysis even in 3D, when B atoms are modeled as fixed 
point-like scatterers, as is known for light waves |23ll24l |. 

It is the 3D version of this scheme that we analyze in 
this paper. After the detailed presentation of our model 
and its practical implementation given in Sec^I we show 



The scatterers are a set of identical particles, whose 
chemical species and quantum numbers will be indexed 
by the letter B, randomly occupying (with filling factor 
p < 1) the sites of a 3D cubic optical lattice. The poten- 
tial seen by the B atoms is produced by a superposition 
of three laser standing waves of common intensity and 
mutually orthogonal linear polarizations along the x, y 
and z axes, 

l/^(r) - Vi'lsiu'ikLx) + sin2(fciy) + sin2(fciz)], (2) 

where Vq^ > is the modulation depth of the lat- 
tice and k^ = l-KJXh is the laser wavenumber. We 
shall denote the lattice spatial period along each axis by 
d = \l/2 = 7r/fci. Multiple occupation of a lattice well 
is assumed to be absent, by use of polarized fermions, or 
by creation of vacancies in a unit occupancy Mott phase 
state U, or simply by choosing p sufficiently small to 
make it statistically irrelevant. 

We choose the lattice depth V^ to be much larger than 
the recoil energy E^ = f{^k\/2mB of the B atoms so 
that the tunneling time of B atoms from one lattice site 
to another (itunnd ~ 1-6 x Wh/E^ for V^ = 50-B^) is 
negligible over the duration of the experiment and the 
disordered spatial pattern of B atoms is static [23 . 

We also require that each B atom is prepared in the 
vibrational ground state of the local lattice micro-trap, 
which can be achieved in practice by raising adiabatically 
the optical lattice on top of a condensate cloud of atoms 
B I2I, or b y app lying Raman laser cooling sideband tech- 
niques [2^ [23 to an optical molasses. This condition 
is crucial to ensure that each A — B scattering event is 
elastic when the A atoms have sufficiently low energy: in- 
deed, energy conservation guarantees that the B atom is 
left in the vibrational ground state after scattering with 
a A atom of momentum k if 



2l2 



2mA 



<.h^^2{V^BE^f/^ 



(3) 



where Lo is the oscillation frequency of a i? atom in a 
micro-trap. 

A last point is to ensure that spontaneous emission pro- 
cesses are negligible for the B atoms. In order to achieve 
large values of V^ with negligible heating of the trapped 
scatterers, we require the lattice to be blue-detuned with 
respect to the strongest transition of the B atoms (in 
blue-detuned lattices, particles are trapped in the min- 
ima of intensity of the stationary light field). Including 
the Lamb-Dicke type reduction factor coming from the 
trapping of B atoms close to the nodes of the laser field, 
one gets the fluorescence rate 



V^ 



3fci 



llJj^ — ujb 2niB0J 



(4) 



where ujl ~ ^b is the atom-laser detuning and Tb is the 
spontaneous emission rate of B atoms. 

For the bosonic *^Rb isotope of rubidium {Xb,D2 = 
780nm, Xb.di ~ 794. 8nm) and an optical lattice tuned at 
Ai = 779nm, only at Inm to the blue of the strongest ru- 
bidium transition at 780nm, at the required lattice inten- 
sity V^ = 50E^ the tunnehng time is itmmoi « 0.7s pTJ 
and the fluorescence rate is Fg^^ ~ 10~*i?^/fi. ~ 3s~^, 
allowing experimental times up to 300ms. The same cal- 
culation, taking for B the fermionic isotope ^"K of potas- 
sium {Xb,D2 — 766. 5nm, Xb.di — 769. 9nm) and an opti- 
cal wavelength A = 765. 5nm leads to ftunnci ~ 0.32s and 

'^ fluo ^S ■ 



B. Model Hamiltonian for the matter wave and its 
limitations 

The matter wave to be strongly-localized is made of 
atoms of another species, that we will label by A. We 
shall ignore interaction effects among these A atoms. One 
way to fulfill this condition in a real experiment would 
be to take spin polarized fermionic atoms: s-wave in- 
teractions are prohibited by the exclusion principle and 
p-wave interactions are very weak at low energies in the 
absence of a p-wave resonance. 

The A atoms experience interactions with the trapped 
species B. At low incoming kinetic energy of a A particle, 
we model these interactions by static contact potentials, 
corresponding to infinitely-massive point-like scatterers, 
each located at the center of a micro-well occupied by a 
B atom: 



N 



V = 2^ S{rA - r^O^i 



rriA 



r^-r,\{\rA~rj\ 



(5) 



i=i 



where the sum is taken over the N scatterers. The effec- 
tive scattering length aoff of a A atom on a trapped B 
atom, when expressed in units of the harmonic-oscillator 
length Oho — \/nJrnBUj, is a function of the dimensionless 
ratios rriB/rnA and a/oho, a being the A — B scattering 
length in free space. The value of a^s and the validity 



condition of our model potential, Eq. (jSj), will be given 
inSec.lTVI 

The A atoms also experience the optical lattice poten- 
tial, with the same spatial dependence as in Eq. (J2J) but 
with a different modulation amplitude Vq. We require 
the optical lattice to be much closer to resonance with B 
atoms than with A atoms, jw^ — lua\ 2> {lol — i^b\, such 
that Vff^ will be much smaller than Vjj^. In particular, 
we impose that 



W' 



«^'-S 



(6) 



SO that, in the absence of B atoms, the A atoms can be 
safely considered as free. In this respect, a particularly 
promising combination is given by fermionic ^Li for the 
species A {Xa ~ 671nm) and ^^Rb for the species B: 
taking Xl = 779nm and a laser intensity such that Vg^ = 
50E^, one finds V^f" = -OME^. If one takes for B 
the fermionic ""^K with Xj^ = 765. 5nm, one finds V(^ = 

-o.oQE^ m. 

We shall therefore neglect the effect of the optical lat- 
tice on the A atoms and take as a model Hamiltonian for 
the matter wave: 



n = nQ + V with Ho 



2mA 



-Ar 



(7) 



We note in passing that in the original Anderson model 
the A particles were instead assumed to be in the tight- 
binding regime, so that the strong localization reported 
in this paper is not stricto sensu Anderson localization. 

To end this section, we briefly discuss two effects not 
included in our model Hamiltonian that may impose lim- 
itations in a real experiment. As we shall see, the produc- 
tion of localized states with a short localization length (of 
the order of the lattice spacing d) requires the use of a 
large and positive value of UeS ~ d, obtained hy a. A — B 
Feshbach resonance. As a consequence, the matter wave 
A has a weakly bound state with a trapped atom B, of 
spatial extension ^ a^g. A first undesired effect is there- 
fore the formation of such A — B dimers. 

A first stage that may lead to a dimer production is 
during the Feshbach ramp of flos from '--^ to ^ d. This 
may be avoided by using a ramping time longer than the 
inverse of the dimer binding frequency, 2m^a^jj/fi, ~ 30/iS 
for our previous example with lithium and rubidium. 
Once flcff is set to ~ d, one may fear that three-body 
collisions A + A + B lead to the formation of a dimer. 
For our model Hamiltonian Eq. (0) the A particles are an 
ideal gas and the dimer formation does not happen: the 
trapped B particle, being replaced by a fixed scatterer 
with no degree of freedom, cannot mediate a ^ — A in- 
teraction. In the opposite limit where the B scatterer 
is supposed to move freely, the rate of dimer forma- 
tion per B atom is jdm = Cdimikpd)^ E^ /h, where we 
used Eqs. (11,12) of [23, taking a dimer binding energy 
e = fi,^/2/ia^g with /i the A — B reduced mass, aoff = rf, 
and assuming that the A atoms are degenerate fermions 



of Fermi momentum kp with kpd < 1; the constant Cdim 
is 6 X 10^^ for the A = ^Li, B = ^'^Rb case, resulting for 
kpd < 1/2 in a dimer formation rate much smaller than 
e.g. the B fluorescence rate T^^^ due to the lattice |30l |. 
The calculation of the actual dimer formation rate in our 
model, taking into consideration the trapping of the B 
atoms, requires the solution of a three-body problem with 
no center-of-mass separability, which is beyond the scope 
of the present work. 

A second undesired effect is the one of gravity. If the 
lattice is arranged to be stationary in a free-falling frame, 
this frame having initially an upward velocity component 
V in the lab frame and finally a downward velocity com- 
ponent V, the overall vertical motion of the lattice in 
such a fountain-like experiment is less than 3 cm for a 
total time of 150 ms. Longer times may be obtained if 
one compensates gravity, e.g. by a using the inflexion 
point of the optical potential produced by a far-detuned 
Gaussian laser beam, or by using electro-optical poten- 
tials 1^ . A drastic solution is of course to perform the 
experiment in a micro-gravity environment 



III. LOCALIZED STATES 

The loffe-Regel criterion Eq. (^ is considered as a nec- 
essary condition to achieve strong localization [3. It is 
simple to check that the setup that we consider can sat- 
isfy this criterion for experimentally reasonable parame- 
ters. As we will see in Sec. IIVI the effective scattering 
length presents confinement induced resonances that al- 
low one to reach the unitary regime for the interaction of 
the matter wave with a trapped B atom, with the maxi- 
mal cross-section a — 47r/fc^. In this case, the loffe-Regel 
criterion reads: 



kd < (Airp) 



1/3 



(8) 



which for a filling factor p = 0.1 yields k < 1/d ~ fci/Tr, 
achievable with sub- recoil laser cooling techniques, or 
with a low density degenerate Fermi gas with a Fermi 
wavevector kp < /cl/tt. 

Since the loffe-Rcgel criterion is not proved to be suffi- 
cient, we numerically investigate in this section the pos- 
sibility for the disordered model Hamiltonian Eq. iQ to 
lead to matter wave localization. In particular, we shall 
find that the unitary regime Oeff = 00 is not the most 
favorable one. 



A. How to find localized wavefunctions ? 

A criterion of strong localization presented by Kramer 
and MacKinnon Q for electrons in a solid consists in 
showing that, at the Fermi energy E = Ep, off-diagonal 
elements of the resolvent G = {E + iO+ — H)^^ in real 
space decrease exponentially with the distance between 
the two considered points. 



For our model Hamiltonian, the calculation of the ma- 
trix elements of the resolvent is straightforward to im- 
plement numerically, using a technique well known for 
scalar light waves in a gas of scatterers [23, UM ■ These 
matrix elements are indeed given in presence of N point- 
like scatterers by 

(rl^lO = 5o(r-r') + ^^x 

X ^,9o(r-r,)[A/-^]^,3o(r/-r'). (9) 

Here go is the propagator in free space of a particle of 
positive energy E = h^k'^ /2mA, fc > 0: 

ik\r—r' 



g,{r-v')^{v\g,\v')^- 



TTiA e 



(10) 



2Trn? |r-r'| ' 
and we have introduced the N x N matrix M: 

]\I= — +M°°, (11) 

flcff 

where / is the identity matrix and Af°° a complex sym- 
metric (not hermitian) matrix with elements defined by 



M: 



ji 



eyi-p{ik\rj-ri\)/\rj -ri\ ii j ^ I, 
ik if j ~ I. 



(12) 



The exact calculation of the resolvent in coordinates 
space is in this way reduced to the inversion of the N x N 
matrix M |3ll |. 

We have implemented the criterion by Kramer and 
MacKinnon for a variable energy E, and indeed we found 
an exponential decay of |(r |^| r') | for sufficiently low en- 
ergies. However, as shown in the appendix ^ this rapid 
decay is not a proof of localization but may be due to the 
fact that £■ is in a spectral gap of the system [33 . 

The most direct way to prove localization is to exhibit 
stationary states that are 'localized' inside the disordered 
potential, that is with a wave function strongly peaked 
inside the scattering medium, decreasing exponentially 
towards the borders of the scattering medium. To this 
end, we use the fact that the wave function 



^{r;ro) = Im{r\g{E + iQ+)\ro) 



(13) 



when not identically zero, is an exact eigenstate of 7i with 
energy E, whatever the arbitrary location of its center 
To 133. Here we take E > so that 0(r;ro) belongs 
to the continuum of the energy spectrum of Ti, like the 
scattering states. A useful expression of is then 



(r;ro) = Aim 



Akli 


■-rn| 


N 


^ik\i 


-•-,■1 


\r- 


-i"o| 


j=i 


|r- 


Tjl 



(14) 



where A ~ ~mA/{2'Kfi^) is a constant factor and dj = 
^Y.M'^^^)3i exp(ifc|r; - ro|)/|r/ - ro|. In practice, we 
choose rg inside the scattering medium; to see if cf) is 
localized or not, one just has to compare the values of 
inside and outside the scattering medium; one can even 
watch how the modulus of decays inside the medium. 



B. Application of tiie proposed technique in ID: 
analytical results 

We test the proposed technique in ID, where analyt- 
ical results can be obtained. We assume that a quan- 
tum particle of mass m interacts in ID with A'' point- 
like scatterers that arc randomly distributed in the in- 
terval X G [—L/2,L/2]. Using the transfer matrix for- 
malism, as detailed in the appendix^ we can obtain ex- 
act expressions for the matrix elements of the resolvent 
■0(a;) = {x\Q{E + iO''')|0) as a function of x. Assuming 
that there is no scatterer in x = 0, we obtain 



^(0) 



m (l + r_)(l+r+) 



(15) 



ih?k 1 — r-r+ 



where t^ is the transmission coefficient of the set of scat- 
terers in the half-space x > oriented from to -|-c». A 
similar expression holds for x < —L/2, see appendix IbI 

Using Furstcnbcrg theorem |34l |. we know that in the 
thermodynamic limit L — > -|-oo, with a fixed density of 
scatterers, the modulus of the transmission coefficients 
tends exponentially to zero, 



U 



\t^ 



\t\ ex e' 



-i/2« 



(17) 



where ^ is the localization length, so that the modulus of 
r_|_ and r_ tends to one. We then find that the wavefunc- 
tion 0(x) = hniplx) outside the scattering medium is of 
the order of \t\. To calculate (f> inside the medium, we 
first take, to zcroth order in |i|, |r+| = |r_| = 1: Eq. p5|l 
then leads to a purely real ip{0) (of order unity), that is to 
a vanishing 0(0). Going to next order, |r±| ~ 1 — |t±P/2 
leads to |(/'(0)| ~ |tp. For a generic value of k, the 
wavefunction (/'(x) has therefore the behavior depicted 
in Fig[2li. 

However, for specific values feg of fc, the above reason- 
ing is incorrect. Assume that k is such that r+r^ is a real 
and positive number. Then the denominator 1 — r_|_r_ of 
Eas. H15ll6|) is extremely small, of order |ip. This leads 
to \(f>{x)\ decreasing exponentially from x = outwards, 
from a value ~ l/|ip to a value ~ l/|i| [33, as depicted 
in Fig[2h- In this case, 4>ix) corresponds to a 'localized' 
state inside the medium. 

Strictly speaking, since we consider disorder over a 
finite region, 4'{x) is not localized: the exponential 
decrease of the envelope stops outside the scattering 
medium, so that 0(a;) is not square integrable. The state 



(t>{x) rather corresponds to a resonance, whose lifetime is 
of the order of h over the energy width of the resonance. 
In practice, can be considered as localized if the res- 
onance lifetime is much longer than the duration of the 
experiment. Assuming that the phase of r_|_r_ varies lin- 
early with k close to fco, one finds that this 'localized' 
state is present on a narrow interval in k of width ex |tp, 
so that its energy width, or equivalently its inverse life- 
time, scales as 



Fcx 



+ \t-\ 



|tp cxe 



-L/i 



(18) 



This illustrates how the proposed technique, defined in 
Eq. ifT^ . gives access to localized wavefunctions. 



where fc > is such that E = h'^k'^/2m, r_ is the complex -^ 



refiection coefficient for the set of scatterers in the half- 
space a; < oriented from a; = Otox = — oo, and r_|_ is the 
complex reflection coefficient for the set of scatterers in 
the half-space x > oriented from x = to a; = -l-oo. In 
the half-space a; > L/2 we also have a simple expression: 



u 


~ltl 




' 








~ltl 


f > 


\ 










bX) 






y^ 


o 
























■^^ 




^^<!}>^ 




^ 
^ 






(a) 








1 



-0.5 




x/L 



0.5 



(16) i 



^ 
^ 



~i/iti y 


2 

~i/itr 


(b) 

\ ~l/ltl 




1 



-0.5 




x/L 



0.5 



FIG. 2: Schematic view of the behavior of the wavefunction 
4>{x) = lin{x\Q\0) in presence of a ID scattering medium of 
length L much larger than the localization length ^, so that 
the modulus of the transmission coefficient for each half of 
the medium jt+|,|i_| ~ |f| oc e~'"^^^ <^ 1. (a) For a generic 
positive energy: the wavefunction decreases exponentially in- 
side the medium, being of modulus ~ |f| out of the medium 
and ~ |ip in the center a:: = of the medium, (b) For specific 
values of the energy, the wavefunction is 'localized' inside the 
medium: its modulus decreases from ~ l/|ip in x — to 
~ l/\t\ for \x\ > L/2. 



Application of the proposed technique in 3D: 
numerical results 



We now present numerical results obtained for a single 
realization of a random potential obtained by a Monte 
Carlo generation of the positions of a finite number N 
of scatterers at the nodes of a cubic lattice with a given 
filling factor p. 

In Fig. O we plot the ratio of the square of the ampli- 
tudes of (/) inside and outside the scattering medium, for 
N ~ 900 scatterers and a filling factor p = 0.1. In or- 
der to avoid a choice that might pick a node of 0(r), we 
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FIG. 3: (Color online) Ratio of the values of 0^ in the center 
and outside the scattering medium (see text) as a function of 
k = {2mAEy /h, with acff — d, for a given realization of 
the disorder with A'' ~ 900 scatterers that occupy the nodes 
of a cubic lattice with 21 sites per side, with a filling fac- 
tor p = 0.1. The value of (jjout is calculated as explained in 
[3bI (red solid curve) or with the extrapolation from the far 
field behavior, Eq. (I20II with r — R — 20d (green dashed 
curve). The energy intervals where the matter wave signif- 
icantly penetrates the scattering medium correspond to the 
narrow peaks in this figure. We have checked that the wave 
function is actually spatially localized in such an energy inter- 
val. The vertical lines mark the locations of the resonances 
obtained by the spectral method: dashed lines for the very 
long lived resonances (F < 6 x 10~^E:^/h) and dotted lines 
for the broader resonances. The temporal decay rate F is 
obtained from Eq. 12311 . 



plot 't'fn/'Pout^ where both numerator and denominator 
are averaged over a few points JBq . The graph reveals 
that the phenomenology is similar to the ID case: one 
has generically (/)?j ^ '/'outJ except for narrow energy in- 
tervals, corresponding to the peaks in the figure, where 
the matter wave can significantly penetrate the scattering 
medium. We have verified for a large number of peaks 
that the wavefunction (f> is indeed 'localized' inside the 
medium. 



We illustrate this phenomenology for a generic value 
of k, and for one that corresponds to a peak in the func- 
tion 4'f^/(f>out taking now a larger number of scatterers 
TV ~ 3400. As can be seen in Fig. 01 in the generic case 
(/>(r) decays essentially exponentially when entering the 
medium, whereas for the specific value of A; the wavefunc- 
tion decays essentially exponentially from the medium 
center towards the outside. What is the associated local- 
ization length ^ ? If one takes ex ex p(— |r — ro|/^), one 
obtains from a fit the estimate ^ ~ d |37| . 




FIG. 4: Plot of 0(r)^ along the straight line passing through 
ro and parallel to z axis, for two values of k, (a) a generic value 
k = 0.3/d and (b) a value k = 0.350134274724/d correspond- 
ing to a peak for (^in/<?!>out as a function of k (in the spirit of 
Fig|3J . Note the similarity with the ID case sketched in Fig[5| 
In (a) the position ro is close to the center of the scattering 
medium: ro = (rf/2, —d/2, —d/2); in (b) it is close to the 'cen- 
ter' of the locahzed state: ro = (15d/2, -3d/2, -lld/2). The 
effective scattering length is given by d/fleff = 1.20530122302. 
To get a clear evidence of the exponential decay of (p, we used 
a larger scattering medium than in Fig.|3] A'^ ~ 3400 atoms on 
the lattice within a sphere of radius 20d, with an occupation 
probability p = 0.1. 



D. Comparison with a spectral technique 

In this subsection, we adopt a different point of view on 
localization jSg , based on the properties of the analytic 
continuation of the resolvent G{z) = l/{z — Ti) from the 
half-plane Imz > to the half-plane Imz < 0, across the 
branch cut Im(z) = 0, Re (z) > corresponding to the 
continuous spectrum of the Hamiltonian. This analytic 
continuation can have poles in z = Eq — ihT/2, where Eq 
and r are positive (see Fig. EJ: these poles correspond 
to resonances. The lifetime of the state associated to the 
resonance is given by F"^. 

The key property of localized states is that they cor- 
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FIG. 5: (Color online) Analytic properties of the resolvent 
G{z) = l/{z — 7i) in the complex plane z. The resolvent has 
poles on the negative side of the real axis, corresponding to 
bound states, and a branch cut on the positive side of the real 
axis, corresponding to the continuum. The analytic contin- 
uation of the resolvent across the branch cut from Im z > 
to Im 2 < may present a discrete set of poles in the fourth 
quadrant: the associated states are the resonances of the sys- 
tem. 



respond to resonances with a decay rate F that tends to 
zero in the limit of an infinite extension of the disorder, 
so that they become in this hmit square integrable sta- 
tionary states of the Hamiltonian. We expect that the 
poles associated to these narrow resonances will leave sig- 
natures on the real axis, in the form of eigenvalues of 
the matrix M with vanishing real part and tiny imagi- 
nary part, for an energy E ~ h?k^ /2mA close to the real 
part Eq of the poles. As opposed to the poles of Q, the 
eigenvalues of M can be calculated in a straightforward 
manner, and in Fig.|31we show with vertical lines the val- 
ues of k for which the matrix M has a purely imaginary 
eigenvalue: the ones with the smallest imaginary parts 
(dashed vertical lines) are in good correspondence with 
the narrow peaks in 0fn/</>out- 

Now, for a value of k within the width of a peak in 
FigO we give an analytical argument relating the spa- 
tial decay of 0(r) to the presence of a tiny eigenvalue mg 
of M of modulus much smaller than the other eigenval- 
ues. In the large r limit, Eq. (|14|l reduces to the far field 
expression 



green dashed curve in Fig|3| We now need to estimate 
the wavefunction inside the medium. Close to the scat- 
tering center located in Tj , one finds cj) ^ ^ Im d^ / 1 r — r^ | . 
This suggests that the dj are also localized, in the sense 
that they decrease rapidly from a center in the medium 
outwards: this we have checked numerically. Then, av- 
eraging spatially over a small volume Z'^, I being smaller 
than the mean distance between scatterers, one finds that 
0?j is at most of the order of {A/l)'^maxj\lm. dj\'^, a value 
reached when Tq is close to the 'center' of the localized 
wavefunction. We reach the estimate 
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maxj|Imdjp 
i* ■ {Im M/k)d 



(21) 



In the right hand side, the first factor has a geometrical 
origin whereas the second one is sensitive to matter wave 
interference effects due to multiple scattering on the B 
atoms: since the matrix elements of Im M/k arc of the 
order of unity, only interference effects can indeed lead to 
a very small expectation value of this matrix. Using ar- 
guments detailed in the appendix[ni we ultimately arrive 
at 



blut 






(22) 



which links the spatial decay of the wavefunction to the 
smallness of an eigenvalue of the matrix M . A useful 
application of this formula is to give the shape of the 
resonances in Fig|2| We linearize the k dependence of 
Too around the value feg such that toq is purely imagi- 
nary: mo{k) ~ f3(k — ko) + ia, where a > 0. Anticipating 
Fig. 13 it appears that the imaginary part of an eigen- 
value of M°°, when tiny for k = kg, remains tiny for 
even lower values oi k, so that the derivative of Immo is 
also tiny; since the real part varies on the contrary over 
an interval of width 1/d, its derivative is not extremely 
small: (3 is essentially real ^^- Eq. (|22|) then leads to 
a Lorentzian shape of the peaks in Fig|31 in agreement 
with the numerics. Since the peaks have a width much 
narrower than /sq, this also leads to a Lorentzian depen- 
dence with the energy E = h'^k'^ /2mA, with a half-width 
at half maximum 1411 : 
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-ikn-TQ 
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-'iA;n-r^ 



(19) 
with n = r/r. As shown in the appendix [Q the sum over 

j in the right hand side is typically (Iml/mJ)^/^ times 

larger than the first term so that the angular average of 

0^ is given, apart from oscillating terms ^ e^''^''/r^, by 



(0(r;ro)^ 



A' 



d * ■ (Im M)d 
2k^ ■ 



(20) 



This allows to estimate 0out by extrapolating this far 
field expression down to r = i?, where i? is a distance 
of order the size of the scattering medium |33], see the 



hT h?kQ ImTO,o(fco) 
2 TOA ReTOQ(fco) 



ex ImTTlg. 



(23) 



When combined to Eq. (|22|) , this leads to a formula nicely 
relating the inverse lifetime F of the localized state to its 
spatial decay: 
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F(x-^(fc = fco)- 



(24) 



If one then assumes a state with a wavefunction localized 
close to the center of the medium and decaying exponen- 
tially as e~^^^ away from the center, Eq. H22|l leads to 
ImTOo ex e~^/^, where L is the diameter of the scatter- 
ing medium, and Eq. H24|) leads to a resonance energy 



thus obeying the same scahng as in 



width r a e-^/«, 
ID, see Eq. (O. 

An important practical consequence of the present 
spectral approach is to give at once the range of values 
of Gcs for which one can hope to have localized states. 
For a given value of k, M will have a tiny purely imagi- 
nary eigenvalue if l/coff is opposite to the real part of an 
eigenvalue of M°° with a very small imaginary part, see 
Eq. Hll|l . In Fig. the eigenvalues of M°° are shown 
as points in the complex plane, for various values of 
the incoming wave number k and for two realizations 
of the disorder with widely different densities of scat- 
terers. We observe that, as the incoming wave number 
k decreases below the inverse of the mean distance be- 
tween the scattercrs p^'^ /d, many eigenvalues acquire an 
extremely small imaginary part and accumulate in the 
region Rem°° ^ —2p^/^/d (see the green dashed lines in 
Fig . 0) , similarly to earlier calculations for light waves 

In a matter wave experiment, this suggests to tune a~g 
to a value close to 2p^^^/d (as we have done in Fig|3J): 
this might be achieved in practice by using a Feshbach 
resonance, as shown in the next section. In this way, 
as k decreases below 2p^/^/d, one obtains a very large 
sequence of values of k such that M has a tiny purely 
imaginary eigenvalue, that is one may have access to a 
large number of localized states. This is illustrated in 
FigEI 

Finally, the representation in FigElis also useful to un- 
derstand what happens in the low energy limit, fc ^ 
It shows that, for each eigenvalue m°° of the matrix M, 
ReTO°° and Iuim°°/k have a finite limit |43j : a numer- 
ical inspection reveals that some of the eigenstates of 
M°°{k = 0), having a tiny value of lim/j^o Imm"" /k, are 
localized. Can these localized states in the zero energy 
limit be accessed in a real experiment ? For a given real- 
ization of disorder, this would require that Ucs be tuned 
exactly to one of the corresponding values of — 1/Re77i°°, 
which is unrealistic. 

In the same way, it would be very difficult to adjust 
acff to hit, at a given value of /c, one of the peaks in 
FiglZI Fortunately, in a real experiment, Ceff is fixed and 
a broad interval of k can be populated by the atomic 
wavepackets; then, FigEI shows that the real parts of the 
eigenvalues of M°° with tiny imaginary parts are increas- 
ing functions of A:, so that they have a chance to cross the 
value — l/ooff for some value of k within the experimen- 
tally populated interval, and to lead to a peak in Fig|2| 
and to an observable localized state in the experiment. 

To conclude this section, it is useful to discuss the num- 
ber of localized states that can be supported by a finite 
size scattering medium, to have an idea of the number 
of atoms that may populate these localized states in a 
real experiment, and to have an estimate of the number 
of states that can be considered as localized for a given 
duration of the experiment. To this end, we present in 
Fig|Hl an histogram giving the number of purely imagi- 
nary eigenvalues of AI per class of inverse lifetimes F. It 
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FIG. 6: (Color online) Representation in the complex plane 
of the eigenvalues m°° of M°° for different values of the wave 
number k of the matter wave. The eigenvalues of M are 
then simply deduced from these eigenvalues by a shift of 
1/aefi along the real axis. Left column: cubic lattice with 
21 sites/side, p = 0.1 and 872 scatterers; the inverse mean 
distance between scatterers is p^'^ /d ~ 0.46/d (same realiza- 
tion of disorder as in Fig.|3Jl. Right column: cubic lattice with 
209 sites/side, p = 10^* and 878 scatterers; the inverse mean 
distance between scatterers is p^'^/d ~ 0.046/d. The green 
dashed lines mark the values Re(m°°) = —2p^'^/d. The real 
axis is in units of p^'^/d, and common values of kd/p^'^ are 
taken in both columns, so as to reveal a possible universality 
in the low p limit. The imaginary axis is in units of A:, as 
justified by Eq. ^. 
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FIG. 7: (Color online) For a fixed positive energy, correspond- 
ing to A; = 0.35/d, ratio of the values of (^^ in the center and 
outside the scatterin g rn edium as a function of the effective 
scattering length a^s ISy] . The dashed line gives the result for 
Oefi = 0. The same realization of disorder is used as in Figf^ 
The inset is a magnification of the region Ocff ~ d. 



is apparent on this figure how an increase of the volume 
of the scattering medium (here by a factor ^ 3.3) leads 
to both an increase in the total number of localized states 
(for a given lifetime) and to the appearance of a tail of 
states of significantly longer lifetimes (here by about two 
orders of magnitude) . 



IV. EFFECTIVE INTERACTION 

We will solve here the two-body problem of a free atom 
A scattering on a single trapped particle B prepared in 
the ground state of a harmonic oscillator potential. In 
the limit of vanishing incoming kinetic energy of the A 
particle, this gives access to the effective scattering length 
introduced in Eq. ^ , which is the key parameter of our 
model Hamiltonian with disorder. 



FIG. 8: (Color online) For a single realization of disorder, 
histogram giving the number of imaginary eigenvalues of M 
per class of inverse lifetimes F^^. The filling factor of the 
lattice is p = 0.1 and the effective scattering length is acs = 
d. Black dashed histogram: A'^ = 872 scatterers within a 
cube with 21 sites per side. Red solid line histogram: N = 
2985 scatterers within a cube with 31 sites per side. These 
histograms were constructed by a dichotomy search of the 
values fco of the momentum fc £ [0, 1/d] such that the matrix 
Af has a purely imaginary eigenvalue; the associated lifetime 
was then calculated with Eq. 1231 . Eigenvalues with negative 
values of F are of course not included (see |4ll['l. The decay 
rate F is given in units of the recoil angular frequency E^ /h = 
hk\/2mA of the species A. In practice, only the eigenvalues 
with small enough values of F are expected to produce an 
observable resonance in (/>fn/<^out as a function of k, see Fig|3 
If the matter wave is made of atoms of ®Li, a duration of 
the experiment of 0.3s corresponds to a minimal observable 
inverse lifetime of F ~ 1Q~^ E^ /fi for an optical lattice laser 
wavelength of 779 nm. 



which represents a spherical matter wave {A) of wave vec- 
tor k and a particle {B) in the ground state of the har- 
monic potential, (l)o{rB) = exp{~r%/2alJ/ (y^aho) , 
with the harmonic oscillator length Oho = (fi/m-sw)^'^. 
The state tJjq has an energy 



A. Calculation of the A-B scattering amplitude 
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H-'k 
2m. 



-huj. 



(27) 



The Hamiltonian for a free A particle and a harmoni- 
cally trapped B one is, in the absence of A-B interaction. 



Hn = — 



2mA 2771b 
It admits the s-wave solution H4 



-mBOj'^rl. 



sm(fcrA) , , . 

ilJo[rA,rB) = — ; 0o(rB), 

kvA 



(25) 



(26) 



We model the A-B interaction with a regularized contact 
potential, which leads to the full Hamiltonian: 



H = Ho + g5{YA-VB) 



d 



d\vA -rsl 



Ir^ -rsl 



(28) 



Here, the derivative is taken for a fixed value of the cen- 
ter of mass position R of the two particles, the coupling 
constant g is expressed in terms of the reduced mass 
/i = TTiAmB/iTTT-A + thb) and the free-space scattering 
length a (relative to the A-B interaction in the absence 
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of the trapping potential) by 



9 



l-i- 



(29) 



The Schrodinger equation Hip = E^p can be reformulated 
equivalently in the integral form 

■>p{rA,rB) = V'o(i"A,rs)+,9 / dpG'B(rA, r^; p, p)V'rcg(p) 

in terms of the two-particle retarded Green's function for 
the non-interacting Hamiltonian 



Ge — 



1 



E + iO+ - Ha 



(31) 



and the regularized part of the two-particle wave 
function, 



V'reg(R) 



d 



d\vA -rsl 



]rA -rB\ip{vA,VB)] 



(32) 



where R = {niAVA + '>^bT^b)/(.'^a + ^^3) and if^^^gi'^) ~ 
Tp-ccgiR) since we consider s-wave scattering only. Insert- 
ing Eq. H30() into the definition of V'reg, we find that the 
regularized part of the wave function satisfies an equation 
of the form 



l/'rc 



I-gO 



aV'o, 



(33) 



where ^^(i?) = Rt}j{R), Vo(-R) = R(t)o{R)s\VL{kR)/{kR), 
and O is an integral operator independent of the scat- 
tering length a. The detailed derivation of this equation 
and the explicit form of O are presented in Appendix IeI 
As we now show, the knowledge of f/'rcg directly leads to 
the value of the scattering amplitude /^ for the A wave. 
Expanding Ge on the basis of eigenstates |k^, n) of Hq, 
and projecting in position space, one is able to calculate 
the integral over the wavcvector k^ so that one is left 
with a sum over the vibrational states of the B particle: 



GB(rA,ri3;p,p) 
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Ak\rA 
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47r \VA - p 
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nirB)4'niP) 



0o(i-b)0o(p) 

(34) 



The low-energy assumption (PJ ensures that the only 
open exit channel corresponds to a B particle in the 
ground vibrational state, and therefore that 



hV 



2mA 2m,A 



■ huj {rix + Uy + 7^2) < 0, V n 7^ 0. (35) 



As a consequence, the terms of Eq. (|34|) involving excited 
states of the harmonic oscillator give exponentially van- 
ishing contributions and can be neglected when A is at 
a distance 3> flho from the center of the harmonic well. 



Expanding |r^ — p\ ~ ta ~ p ■ va/ta and substituting 
into Eq. l|Sn|) . we get 



V'(rA,rB) 



sin(fcrA) J. exp(ifcrA) 
+ Jk 



Mrs), 

(36) 



krA ' TA 

where we have introduced the scattering amplitude 

Ik = -a / dp — j^ (/)o(p)V'rcg(p)- (37) 

In the limit of zero-energy, the scattering amplitude de- 
fines the effective scattering length Ceff through 



_ ,. „ mA 
Ooff = - hm jk = a 



dp0o(p)Cg (P)- (38) 



We have solved numerically Eq. H33(l . and in Figs. El- 
111! we plot Ooff as a function of l/a for different values 
of m^s/naA- We choose m^B/mA = 0.15,1,6.67, corre- 
sponding to the physical cases of a mixture oi A = *°K 
and B = ^Li, a mixture of two different internal states of 
atoms of the same species, and a mixture of ^4 = ^Li and 
B = *°K, respectively. It is apparent that a^s presents 
a series of intriguing resonances, which can be used to 
tune Ooff to a high value and whose physical original are 
discussed in illVBI 

In Figs. I lOl and 1 121 we also plot, in the case ms = tua, 
the behavior of the effective range re . It is defined as a 
coefficient in the low k expansion of the inverse scattering 
amplitude: 



/. 
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ik 



Jp-12^ 



(39) 



The replacement of a trapped B particle by a fixed point- 
like scatterer, as done in Eq. ©, is allowed when the fc^ 
term in the above expansion is negligible, i.e. when 



kelfc^ < la J 



-ik\ 



(40) 



In the ideal regime for matter wave localization, fc < 
a^Jjj ^ 1/d. Since r^, is generally of order aho ^ c^, this 
condition Eq. (|40|l is satisfied. Note that rg diverges when 
Q-off — ^ 0, a generic phenomenon: expanding fk rather 
than /^ in powers of fc. 
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(41) 



one generically expects that ^ga^ff has a finite limit when 
Ocff -^ 0. This was demonstrated analytically both for a 
square well |45j and a van der Waals interaction potential 
[iq . In this limit, Eq. H40|) is thus violated for any finite 
A:; however, this is not an ideal regime to obtain matter 
wave localization. 



B. Resonances of the effective scattering length 

We now show that the effective coupling constant ex- 
periences an infinite set of resonances due to the presence 
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FIG. 9: (Color online) Effective scattering length flofr (con- 
tinuous line) as a function of a~^ for mB/rriA ~ 0.15 (trapped 
Li and free K). The vertical lines (red) mark the positions 
of the resonances (a step of Oho/a = 0.01 is used to sample 
the curves, and some of the resonances are too narrow to be 
seen on the graph). 
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FIG. 10: (Color online) Same as Fig. |n|for two particles of 
equal mass. The blue dashed line is the effective range Ve, and 
the magenta dotted line is reborn ~ —{^/mA)a\a/<i, obtained 
in the Born approximation by replacing i/;rcg(p) with tpo{p, p) 
in Eq. 1371 . A green dot at a = — l.Gaho marks the position 
where a^s ~ d \i Vo^ = 50£f . 



FIG. 11: (Color online) Same as Fig. I^lfor mB/mA ~ 6.67 
(trapped K and free Li). The blue dashed line is the 
Born approximation aoff^Bom = arriA/ p-, obtained by replac- 
ing ^r^=°(p) by 4,o{p) in Eq. (EHJ- 
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FIG. 12: (Color online) Energy dependence of the scattering 
amplitude for two particles of equal mass at a = —l.Gaho 
(green marker in Fig. IIOII . The linear fit at low energy, 
Re{l/fk) = — l/ffleff + refc^/2, yields a^s = 8.3aho and 
fe = 1.135aho (units of a^^ and a^^ on the horizontal and 
vertical axis, respectively). 



|e,) (see Eas. HE18IE19p ). Wc then rewrite Eq. ^ as 
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-A, 



(42) 



of the external confinement. In the limit fc — > 0, the re- 
solvent has real matrix elements in position space, see 
Eq. (|34|l : so does the symmetric operator O, which then 
admits real eigenvalues A^ and orthonormal eigenvectors 



where the scalar product of functions of the single vari- 
able R is defined as {u\v) = 47r J„ dRu*{R)v{R). This 
means that a singularity in Ooff is expected whenever the 
denominator of Eq. (|42|l vanishes, i.e. whenever the in- 
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verse free space scattering length 1/a equals 27rh'^Xi/fi, 
provided that the non-interacting wave function has non- 
zero overlap with the corresponding eigenstate |ei). 

Confinement-induced resonances have been been ana- 
lyzed theoretically in different contexts, such as ID wave 
guides [43, liMligl , 3D optical lattices [50| , and quasi-2D 
condensates |51|- In ID wave guides the effect is partic- 
ularly remarkable: due to the presence of the transverse 
confinement, a contact potential acquires a bound state 
for any value of the 3D scattering length (while in free 
space the contact potential has a bound state only for 
a > 0) . Confinement-induced modifications of two-body 
scattering properties have very recently been observed 
experimentally by the Zurich group in ID waveguides 
|53 | and in 3D optical lattices |53 • In most of the cited 
papers, the underlying translational symmetry and the 
harmonic nature of the confinement permit the factor- 
ization of the center-of-mass motion: this in turn implies 
that a single confinement-induced resonance can exist, 
since only one state in the closed channel is coupled to the 
open channel [ig- When this factorization is not possi- 
ble, as in our setup or in the case of a ID wave guide with 
anharmonic transverse confinement |49| , an infinite set of 
states in the closed channel has non-zero coupling to the 
open one, and an infinite number of resonances appears. 
In practice however only a few of them may be resolved 
and relevant in an experiment since they become increas- 
ingly sharper as Oho/a becomes larger (i.e. KfAolej)!^ — > 
for large i). The position of these resonances can be pre- 
dicted analytically in various limits, depending on the 
sign of a, as we now discuss. 



1. Position of the resonances for a > 

When a > 0, the pseudopotential admits a bound 
state in which the two particles can "sit" for a variable 
time, forming a molecule that oscillates in the harmonic 
well. To understand this point, we rewrite the two-body 
Hamiltonian 12811 as 



— - + --mBCJ R -— — + g-5{Y)- — {r-) 

2[mA + mB) 2 2/.i or 



1 "^^^ c.V-/.c.^R.r 

2 niA + TUB 



(43) 



and treat the terms in the square parenthesis as a per- 
turbation. The unperturbed part admits the factorized 
eigenstates ^„^;=o(R7r) = (/)„j=o(R-)x('') that describe 
a bound molecule with internal wavefunction x(^) = 
exp(— r/a)/V27rar and centcr-of-mass in an eigenstate 
of the harmonic oscillator of angular momentum I = 
and radial quantum number n > 0. Since both the ini- 
tial state and the Hamiltonian are spherically symmetric, 
conservation of angular momentum allows only Z = in- 
termediate molecular states. Within this unperturbed 
approximation, Oeff diverges each time the energy of the 



oscillating molecule corresponds to the ground state en- 
ergy of the pair of atoms, i.e. at the values of a = Ores 
that satisfy 



2n+ -] huj 



niB ft 3 
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(44) 



As can be seen in the upper part of Fig. 1131 this formula 
describes the position of the resonances with a > in a 
wide region of the graph, since corrections to it arc only 
O(a/aho)' M- 



2. Existence of a resonance for a < 

For a < and for a large enough rriA/mB ratio, the 
presence of at least one resonance can be demonstrated 
by a very simple variational argument performed at the 
unitary limit 1/a = 0. A contact potential characterized 
by the scattering length a can be replaced by the sim- 
ple boundary condition ip{Il, r) = C(R) {r~^ — a~^) + 
0(1), where C is an arbitrary function of the center of 
mass coordinate R. The s-wave Ansatz 



^{R, r) = AAcxp (-i?VA2) ''''^^ '^'^'"'^ 



(45) 



with A and a variational parameters, satisfies the bound- 
ary condition imposed by a unitarity limited contact po- 
tential, i.e. characterized by 1/a = 0. Its energy can be 
calculated from the Hamiltonian (|43|l : the term R • r has 
a vanishing contribution when averaged over this state, 
and the variational energy assumes the minimum value 



£;° 



3, / rriB , 1. I rriA 

-nuj\ h 7: fti^ 1 / ■ 

2 V ruA + ruB 2 V tua + ruB 



(46) 



We might imagine continuously tuning a from 0^, where 
no bound state can exist, towards —00. If E°° < 3fta;/2 
(i.e. UmA/mB > 9/16) abound state for 1/a = is guar- 
anteed by the former Ansatz, and at least one resonance 
for floff must exist in the a < region. 

The following more elaborate calculation allows to 
prove the existence of a resonance for a < 0, for an ar- 
bitrary mass ratio Tn^/ms- One performs a variational 
calculation directly on the integral operator 0'^=°, tak- 
ing as a variational function f{R) ~ Rcxp{—R^/2a^^). 
Using Eas. (|E18IE19|l and performing Gaussian integrals 
we obtain in harmonic oscillator units: 



(/|6«|/) = [f'" 



{[(. 



dr \(1 



-3/2 



+ a-aexp(-r)]~^/H. (47) 



— [T + a ~ a exp 

From the inequality exp(— r) > 1 — r valid for any r > 0, 
we conclude that the integrand is a negative function, so 
that (./|0'''=°|/) < 0. Since the trial wavefunction f{R) 
is proportional to ipQ^'^ (R) = i? 0o (R) , we conclude that 
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O*"'"" admits at least one eigenvector je^) with a non- 
zero overlap with ip^^^ and with a negative eigenvalue 
Ai < 0. The identity Eq. (|42() then implies the existence 
of a resonance in Ucs for a < 0. 



3. Position of the resonances for a < 

When a < and mB/rnA ^ 1, the position of 
the resonances can be found with the aid of the Born- 
Oppenheimer approximation. For a fixed position of the 
massive particle A^ one calculates the energy of the B 
particle, which then constitutes an effective potential for 
the A particle. Restricting for simplicity to the mean 
field regime \a\ <^ aho, one can assume that the parti- 
cle B remains in the ground state of the well, thereby 
creating an effective Gaussian attractive well for the A 
particle. Hence the effective Hamiltonian for A: 



10 F 



Hcs 



2mA 



-Ar. - 



2TTn\a\ expi-rj^/alJ 



TUB 



{^/^:a^o) 



(48) 



As argued above, the characteristic range of the potential 
is of order Oho- This Hamiltonian can be easily solved nu- 
merically, and predicts a divergence of Oeff whenever the 
combination (|a|/aho)('TiA/'7iB) equals the critical value 
for the appearance of a new bound state (see Fig. El 
dashed lines in lower graph). 

In the opposite limit mB/mA 3> 1, there is one res- 
onance left on the a < side. It is intuitive that its 
position flics tends to — oo in this limit, the B particle be- 
ing then perceived by A as a fixed scatterer of scattering 
length a, for which the resonance is obtained for a = — oo. 
At /c = 0; one then expects that iprcg{R) oc (f>o{R). To 
formalize this intuition, we expand the integral operator 
O in powers of the mass ratio a = mA/iTiB —^ 0: from 
Eas. HE18IE19p . we get 



O 



fe=0 



Y^^^do + a^^^di + 



(49) 



It is then possible to check analytically, by calculation of 
Gaussian integrals, that one has exactly Oo|0o) = 0; with 
4>o{R) = R(f>o{R). Using perturbation theory, we obtain 
the series expansion of the lowest eigenvalue of 0'^^° : 



Ao = a5/2(0o|Oi|0o) + O(«'/' 



a5/2 
ttV2 



+0(a^/2) (50) 



in harmonic oscillator units. The resulting lowest order 
expression for the resonance position is 



Oho _ _ /2 / ™^ 
flics yruB 



3/2 



l + O 



niB 



(51) 




and is shown as a dashed line in Fig ll3l |5E 



V"^A 



FIG. 13: (Color online) Position of the broadest resonances 
for positive (upper graph) and negative (lower graph) a as 
a function of niB/mA- The dashed lines are the theoretical 
predictions: for a > they are given by Eq. I|44|l : for a < 
and ma < niA they are given by the Hamiltonian in Eq. 1481 ; 
for a < and tub > tua it is given by Eq. II5H . In the upper 
graph, from top to bottom symbols correspond to n = 4, 3, 2, 1 
in Eq. 1441 . The dotted vertical line indicates the value of the 
mass ratio mB/rnA when A — Li and B = Rb. 



V. EXPERIMENTAL OUTLOOK AND 
CONCLUSIONS 



In this paper, we have shown numerically that the way 
of producing a disordered potential for matter waves pro- 
posed in J23|, i.e. the use of atoms randomly trapped 
at the nodes of an optical lattice, indeed leads to the 
appearance of exponentially localized states in a three- 
dimensional geometry. Our numerical method allows to 
directly compute the wavefunction of the localized states; 
it is based on the fact that the matrix elements of the re- 
solvent of the Hamiltonian can be calculated extremely 
efficiently for the interaction of the matter wave with 
point scatterers, a fact already used with success in the 
context of light locahzation (H, 0, E3- The method 
also allows to obtain analytical results in a straightfor- 
ward manner, such as a check of the existence of localized 
states in ID from Furstenberg theorem or the derivation 
of a link between the inverse lifetime of a localized state 
and its spatial decay in 3D, for a finite spatial extension 
of the scattering medium. 

The main physical result is that numerous long-lived 
localized states appear for a wave number of the mat- 
ter wave smaller than the inverse of the mean separation 
between scatterers, when the effective scattering length 
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Oeff of the matter wave on a trapped atom is positive 
and of the order of the mean separation of the scatterers. 
For the numerical examples of this paper, with a 10% 
occupancy of the lattice sites and for a wavenumber of 
the matter wave as large as ^ 0.5/(i, there are localized 
states with localization lengths that can be as small as the 
lattice period d, usually a sub-micron quantity. This ex- 
tremely strong localization allows to have very long-lived 
localized states (lifetime larger that 10^ inverse recoil an- 
gular frequency of the matter wave) even in disordered 
samples with a radius as small as 10 lattice periods. By 
a full solution of the scattering problem of a free atom 
with a harmonically trapped one, we have shown that 
the required large values of UcS ^ d can be obtained by 
using an inter-species Feshbach resonance, and we have 
characterized intriguing confinement-induced resonances 
that appear in this two-body scattering process. 

How to proceed in a real experiment to get evidence 
of these localized states ? A possibility is to extend to 
matter waves what was proposed for light in [2J|: one 
introduces the matter wave wavepackets inside the scat- 
tering medium at a low value of Ccff, then one tunes Coff 
to the desired high value and one lets the matter wave 
evolve in presence of the scattering medium (but in the 
absence of an external trapping potential). After an ad- 
justable time T, one measures the number of remaining 
matter wave atoms A',om(''') in the scattering medium. 
Since the component of the matter wave wavefunction in 
localized states decays exponentially in time with very 
weak rates F, the function A'rem(''') should have a long 
tail, as compared to the case of a purely ballistic or even 
diffusive expansion [Sg. A further check that this long 
tail has a decay rate varying exponentially with the size 
of the scattering medium would be a very convincing ev- 
idence of strong localization |53| • 



of a periodic system. 

To illustrate this statement for our system, we have 
enclosed our lattice of scatterers in a box of side L, im- 
posing periodic boundary conditions on the walls of the 
box. This amounts to replacing Eq. H1U|) by the particle 
propagator satisfying the correct boundary conditions: 



grir) 



2m y- 

q 



„jqr 



(Al) 



with q = 2TTn/L and n £ Z^ (a triplet of integers). In- 
deed we found that, for the ensemble of scatterers used 
in Fig. El the ground state of the system, once en- 
closed in a box of side L = 23d (slightly larger than 
the scattering medium) , is characterized by a wave num- 
ber /cniin — 0.7202(i~^. The exponential decay shown at 
k = 0.3d~^ by |(r|tj|r')| is therefore simply indicating 
that at such low energy no state can exist deep inside the 
medium. In a scattering experiment, we might imagine 
a plane wave coming from infinity that scatters on the 
trapped B atoms: if fc < fc„ii,i the incoming wave under- 
goes total reflection, and inside the random medium only 
penetrates an evanescent wave, that decays exponentially 
from the boundary of the medium towards its interior. 
This example clearly points out that an exponential de- 
cay of |(r|C7|r')| is not a sufficient criterion to prove 
localization in our system, since it does not guarantee 
the existence of states deep inside the random potential. 
In the criterion of Kramer and MacKinnon, introduced 
for a solid, the energy is taken equal to the Fermi energy, 
with the assumption that the density of states does not 
vanish at the Fermi energy. This criterion then turns 
out not to be practical in our case, since it requires a 
diagonalization of the Hamiltonian in order to calculate 
the density of states. 
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APPENDIX A: SPATIAL DECAY OF 

{r\g(E + iO+)\r') 

As illustrated in Fig. ^1 in our model the off-diagonal 
matrix elements of the resolvent {r\0{E + iO^)\r') decay 
exponentially in |r — r'| at low energy E. Is this a signa- 
ture of localization ? No, because this decay may be due 
to the fact that the energy i? is in a spectral gap of the 
system. And this occurs as well in the forbidden bands 



APPENDIX B: ANALYTICAL RESULTS IN ID 

We consider a ID free space geometry with N Dirac 
scatterers located in distinct positions xi < . . . < xn- A 
quantum particle of mass m interacts with the N scatter- 
ers, with a coupling constant g. Using the transfer matrix 
formalism we calculate the resolvent Q{z), for z = E + ir], 
E real and 7y > 0, which will give access to stationary 
wavefunctions in the limit ?/ — > 0"*" . Taking for simplicity 
xo = different from the positions Xi, we use the fact 
that ipi^) = {x\Q{z)\x[) = 0) solves Schrodinger's equa- 
tion with a source term: 



2m dx^ 



N 

g^S{x-Xi) 



ij{x)=6{x). (Bl) 



Over a position interval containing neither one of the 
Xi nor xq = 0, ipi^) is an eigenstate of d'^/dx'^ so that, 
introducing the unique kc such that z — h^k^/2m with 
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' k=0.3d 
■^^ k=3.0d" 

2 



lgo(i-)l' 
exp(-lr-r(,l/0.24)/100 



1 



Irl 



FIG. 14: (Color online) Decay of the off-diagonal real space 
matrix elements of the resolvent as a function of the distance 
from the center of the cloud of scatterers (in units of the 
lattice spacing d), averaged over 100 different realizations of 
the random potential. The « 4600 scatterers (aoff = 0.3d) are 
distributed in a cubic box of side 21d, occupying each site of 
the lattice with p — 0.5. The y-axis is in units of [h^ /2mA)~'^ ■ 



Imfcc > 0, one has: 

X < Xl 

Xj < X < i'j+l 

Xs < X < 

< X < Xs+l 

Xj' < X < Xji^i 

Xn < X 






ip{x) = A+e''^"'' 
tPIx) = A+e^''-^ 
i^lx) = A+e^''-'' 

iIj{x) = B+e'''-^'-' + B-,e-'^-^ 
V'(x) = B+e"^"'-' + B^e-'^""" 



B: 



-ikcX 
-ikcX 
-ikcX 



where we assumed that the first s scatterer positions 
Xl, . . . ,Xs are < and the N — s other ones x^+i , . . . , xn 
arc > 0, and where j runs from 1 to s — 1 and j' from 
s + 1 to iV — 1. All the unknown coefficients A^ and 
B^ shall now be determined from boundary conditions 
obeyed by ip{x). 

A first boundary condition is that tpi^) should not di- 
verge exponentially in x = ±oo. This imposes A'^ = 
and B'^ = 0. 

The other boundary conditions originate from the fact 
that each Dirac distribution in Eq. IJB1|) introduces a dis- 
continuity of the first order derivative of tp{x), whereas 
4>{x) remains continuous. Integrating Eq. IJB1|I over 
an infinitesimal interval around Xj leads to tjj'{x'^) — 
ip'{x~) = 2mg'ip{xj)/h'^ with ip{x'^) = ^{x~) = ^{xj). 
These two equations allow to connect the unknown co- 
efficients on the right of Xj to the ones on the left by 
a two-by-two transfer matrix. For j = 1, . . . , s wc thus 
have 



4 \-p( 4-1 

4/ ^Uj-i 



(B2) 



and for j = s + 1, . . . , A'': 



' ' BJ-. 



(B3) 



A simple calculation leads to the explicit expression 



■* > 17 e^''^'=^3 1 -I- Z7 



(B4) 



with 7 = mg/ih^kc). We then introduce the two fol- 
lowing matrices, one associated to the first s scatterers, 



A=Ps...Pi (B5) 

and the other one to the N — s last scatterers: 



b = p: 



■ P, 



N 



(B6) 



They allow to express the coefficients Af in terms of Ag 
and the coefficients B^ in terms of B 



N- 



At 
A7 

Bt 
B7 



A 



= B 



B 



N 



(B7) 
(B8) 



The last two unknowns, A^ and i?^, are obtained from 
the boundary conditions imposed by 8{x) in Eq. (|B1|I . 
V'(O-) = V'(0+) and V''(0+) - V''(0") = 2mlh^, which 
imposes two equations on the coefficients Bf , A^ , 



At = Bt- 
A- = B- + 



m 

ih^kc 
m 

ih?kc 



(B9) 
(BIO) 



Combined with Eq. (|B7|1 and Eq. 
following system: 



B8I). this leads to the 



Bn 
B21 



-A12 

^Al2 



Bt 

An 



N - 



ih?kr 



which can be solved explicitly: 



A^ = 



ni 



Bn + B21 



B 



ih^kc B11A22 - B21A12 

m A12 + A22 



N 



ih-^kc B11A22 - B21A12 ■ 



(Bll) 

(B12) 
(B13) 



We then proceed with the limit z tending to a real and 
positive energy E: 77 — > 0+ and kc ^> k = {2mEy^'^ /h. 
The Pj then become physically meaningful transfer ma- 
trices in the SU{1, 1) group. In the main text we intro- 
duced the refiection and transmission coefficients r+, i+ 
of the last N — s first scatterers on the axis oriented from 
X = —00 to X = +CX3. This means that there exists a sta- 
tionary solution of the usual Schrodinger equation equal 
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to e*'^^ + r+e "^^^ for x < Xg+i and equal to i+e*'"'^ for 
X > xjv- By definition of the transfer matrices: 



One has to calculate the vector d to fully determine 
(r; ro), see Eq. (|14|) . This vector solves the linear system 



B 







1 

r+ 



(B14) 



Similarly, introducing the reflection and transmission co- 
efficients r_ , t_ of the first s scatterers on the axis ori- 
ented this time from x = -l-oo to a; = — cx), we imply the 
existence of a solution of Schrodingcr's equation equal 
to e"*'""^ -I- r-e^^^ for x > Xs and equal to t_e^*'^'^ for 
X < xi, which imposes 



A 



1 



(B15) 



This allows a physical interpretation of the coefficients of 
A and B: 



Al2=r^/t^, ^22 = l/t- 

S21 =r+/t+, 6ii = lA+ 



(B16) 
(B17) 



This leads to very simple expressions for the Green's 
function '0(x) out of the scattering medium and over the 
interval between two scatterers containing xq — 0: 

X <Xi: ^{x) = -r^ -r^ ^e ^"^ 



Xs < X <Q : ip{x) 
< X < Xs+i : "0(2^) 



ih?k 1 — rjf-V- 



1 — r+r- 



x>xn: i^{x) = — — -^ J_^+^kx^ 



Md ^ s 



(C2) 



with the source term Sj ~ — exp{ik\rj — ro|)/|rj — ro|. 
As the eigenvalue mo of M is the only one to be extremely 
close to zero, we take the approximate expression 



M" 



1 
mo 



K'oXfSI 



This leads to 



d ~ v"q {vq ■ s) . 



(C3) 



(C4) 



To have access to an estimate of 4>ont j we have to cal- 
culate the expectation value of the imaginary part of the 
matrix M on the vector d. Since M is symmetric, both 
the real part and the imaginary part of M are hermitian 
matrices, with real expectation values, so that 

d* ■ (Im M)d = Im (d* ■ Md) = Im {d* -s). (C5) 

Then using the approximation Eq. (|(M|I we obtain 

'{vh* ■ s*){v'o* ■ s)' 



d* ■ (ImM)d~ Im 



(C6) 



An immediate application of this result is that the quan- 
tity /(n) = J2j rfje"''^"''"^, where n is a unity vector, is 
typically much larger than unity. The calculation of the 
average of |/p over the unit sphere indeed leads to 



{\f? 



d* ■ (ImM)d 



(C7) 



APPENDIX C: AN APPROXIMATE RELATION 
FOR SOME EXPECTATION VALUES OF ImA/ 

This appendix is useful to derive Eq. (|22|) , an equation 
which relates the spatial decay of 0(r; Tq) to the smallness 
of an eigenvalue of M. 

Let us assume that we are at a positive energy E such 
that the matrix M has an eigenvalue mg extremely close 
to zero, much closer to zero anyway than all the other 
eigenvalues. Let vq be the associated eigenvector of M. 
The corresponding adjoint vector is an eigenvector of M^ 
with the eigenvalue ttiq . Since the matrix M is (complex) 
symmetric, M^ is simply M*, the complex conjugate of 
M, so that one may take as adjoint vector the complex 
conjugate vq* of iSq. The imposed normalization condi- 
tion is then 



{Vq\vq) =-(771^ = 1. 



(CI) 



For compactness we use here Dirac's notation, even if M 
does not act in a Hilbcrt space. 



Since |/| ^ 1 in the low niQ limit [63, it is correct to 
neglect the term £-''="■'"0 in Eq. H19|l as was done in the 
main text. 

To have access to an estimate of (j)in, we have to cal- 
culate the maximal value of all the |Imdj|. Let us call n 
the index such that |Ini(i„| is the biggest one. According 
to the approximation Eq. ljC4p . we then have 



Im d„ ~ Im 



V0,ni'IJ0 



niQ 



(C8) 



where wo,n denotes the component n of the vector wg. 

The expressions Eas. (|C6IC8|l greatly simplify when one 
chooses a position Tq that tends towards r„. The fact 
that the results shall not depend on this specific choice 
of ro is established in the appendix^ In this limit, all 
the components Sj of the source term s are negligible as 
compared to s„ ~ — l/|ro— r„| sothatuo-s ^ — wcn/lro — 
r„| and so on. Then Eq. (|^ reduces to 



^lut 



klf 



Im 



'^O.ri 



(C9) 
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As shown in the appendix^ the components of the vec- 
tor w5 Eire real in the hniit of a vanishing JttioI so that Vq „ 
may be pulled out of the imaginary part. Since the \vo,j P 
decrease roughly exponentially at large distances jr^ — r„| 
over a length scale b of the order of the mean scatterer 
separation, as was known from studies of light localiza- 
tion [23, IM E3 j the normalization condition Eq. (jCip 
leads to Wq „ ~ l/{pb^) < 1, where p is the mean scatter- 
ers density. We then get Eq. l|^ . 



APPENDIX D: ON THE FACT THAT SOME 
QUANTITIES ARE ALMOST REAL 

We consider here a value of k such that the matrix 
M has one (and only one) eigenvalue mo of extremely 
small modulus. We then show that the corresponding 
eigenvector vq of M, normalized as in Eq. (|Cip . is close 
to a vector with real components, a property used in 
the Appendix [0 and necessary to obtain the equation 
Eq. ^. 

First we give a physical argument. Starting from the 
approximation Eq. IjCSp . and keeping terms only to lead- 
ing order in l/rrio, we obtain from Eq. (|14[l : 



(r;ro) ~ -Aim 



F(r)F(ro) 



mo 



where 






r - r, 



(Dl) 



(D2) 



Now let us assume that the phase of -F'(ro) (modulo tt) 
changes significantly (that is in a way not tending to 
zero with ttiq) when ro is varied. This means that the r 
dependent wavefunction (j)(r; Tq) spans a subspace of the 
Hilbcrt space of dimension 2 when ro is varied, a subspace 
generated by the wavefunctions Rei^(r) and Imi^(r), or 
equivalently by F{r) and F*{r). This leads to the phys- 
ically suspect situation that two independent localized 
states can be associated to a given resonance. 

To avoid this suspect situation, one is led to the as- 
sumption that the function F(r) has a constant phase 
(modulo tt). Taking the limit r —^ Tj, this implies that 
all the vqj have the same phase (modulo tt). From 
the normalization condition Eq. (|Cip we conclude that 
all the vqj arc real, apart from small terms tending to 
zero with toq- Furthermore, this leads to the conclu- 
sion that the function F(r) has to be real everywhere, 
apart from terms that tend to zero with mo |64| . A con- 
sequence that is important for the appendix [CI is that 
Vq* ■ s ~ vo* ■ s* = — -F(ro)* are approximately real, so 
that the result Eq. (|C9|) docs not depend indeed on the 
specific choice of tq. 

Now we give a mathematical argument to show that 
Vq is almost real. We start from MiJq = moVo and we 
split the matrix M in a real part and an imaginary part, 



M — Mfi + iMj. Since M is symmetric, both M^ and 
Ml are real symmetric. Taking the squared norm of the 
identity MhVq = (mo/ — iMj)vo, where / is the N x N 
identity matrix, leads to 

WMrVoW" = vh* ■ [Mf+Rc{ml)l] vh, (D3) 

where we used the fact that 



v^* ■ MiVq = Im (mo)t;o* • vq, 



(D4) 



that can be proved as in Eq. IJC5|) . Next, we expand vq 
in the orthonormal eigenbasis of Af/, vq = ^^ Calm/.a), 
so that 



Vo ■ 



A//Vo = ^|c„|2m2,^. 



(D5) 



It remains to use the following property of M/, valid for 
an arbitrary vector x: 



X ■ Mjx = k I — — \ yxjC 
3 



E- 



ikn-Vj |2 



(D6) 



to show that < m/ „ < Nk JG^]- As a consequence, 
"^| Q, < Nkmja- Using the expansion Eq. IJD5|) and 
Eq. (jD4|l . one gets an upper bound on iJj* • M'Jv'q so 
that 

llMflvBIl' < [iVHmmo + Re(mg)] \\vhf. (D7) 

In the large scattering medium diameter L -^ +00, 
an exponential decrease of the imaginary part of mo is 
expected, ~ exp(— L/^). Close to the center of a reso- 
nance peak, |Re77io| ~ ImTTio, so that ||MflWo||/||z;o|| = 
O [A^fcexp(— L/2^)]. On the contrary, the density of 
states for the spectrum of M^ is not expected to be expo- 
nentially peaked: the spacing between successive eigen- 
values is expected to scale as 1/A^ at most. This is ap- 
parent in FiglSJ and we have checked it numerically for 
increasing numbers of scattering centers. This implies 
that Vq, when expanded on the eigenvectors of Mr, pop- 
ulates essentially one eigenvector of Mr, the one with 
the eigenvalue nearest to zero. Since this (unit norm) 
eigenvector of Mr is proportional to a real eigenvector, 
we deduce from Eq. IjCip that the components of vq are 
almost real, at the O(Immo)^/^ accuracy level. 



APPENDIX E: INTEGRAL EQUATION FOR Vrog 

To calculate the non-interacting two-particle Green's 
function numerically, we adapt a technique used in [66j. 
The Feynman propagator K associated with the non- 
interacting Hamiltonian Hq |67J, 



Kt{vA,YB;PA^pB) 



VA,VB e-'^"*/'' pj„ps) (El) 
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can be factorizcd as the product of the two factors, Kt 



Kf^Kj^ , the first term describing a free particle, 



K^{^a]Pa) =e 



-i7r3/4 



/j7M_\3/2 f iruA \VA- Pa? \ 

\2TTntJ '^^^ \ 2n t J 

and the second one a particle in a 3D harmonic oscillator, 






mBijJ 



3/2 



of x{t) = 0. In the latter expression, Eq. (|E7|I . the con- 
tribution of the neighborhood of t = diverges as 1/r for 

e*-3/4 J^ (2i)3/2 J^ 

'^M^+X 'H^M^^- 



exp 



V ^ 



27rfi,| sin(a;i)| 



1 /¥ 



ji7r/4 



r V 2 



dt 



^I+Pl 



2 ia,n{ujt) sin(Lji) 



(2t)3/2 



(E3) 



(E9) 



We finally find 



where (j){t) = exp(i7rn/2) for mr < tut < (n+l)7r. Setting 
Kt = for i < 0, Ge is obtained as the Fourier transform 
oiKt, 



Gi 



^^g«(£+«0+)t/ft^ 



(E4) 



For simplicity, we introduce dimcnsionless variables by 
expressing quantities in harmonic oscillator units and, 
even though the derivation has been carried out for a 
generic mass ratio, we restrict ourselves here to the spe- 
cial case ruA = itib = rn (see end of appendix for the 
case k —^ with an arbitrary mass ratio mA/'niB)- 

In order to find the equation satisfied by V'reg, we 
rewrite Eq. 1)30(1 as 

ip{rA,rB) = il)o{vA,VB)+gAcg{R) / dpCs {rA,rB] p, p) 

+ 9 dpGE {rA,rB; p, p) [Acgip) - iprcgiR)] (E5) 

and let the two particles approach each other. We in- 
troduce here the center-of-mass and relative coordinates 
R — {rA + rB)/2 and r = r^ — r^. Let us turn our 
attention to the first integral appearing in Eq. I)E5(I : 



'^=° V-rcgli?) {gF,{R) - ^) + o(l) (ElO) 

that separates out the expected 1/r divergent contribu- 
tion and defines 



1 



^'^^^^j2^Jo "^^ l^|sin(t)+tcos(i)|3/2-(^^^ 

(Ell) 

with Wo = i?2(cos(t) -tsin(t)/2- 1) / (sin(t) -f-<cos(t)). 
The term (2t)~3/2 regularizes Fi in the neighborhood of 
i = 0, and we have taken the limit J — > 0. 

Let us now consider the remaining term appearing 
in Eq. HE5|) : when r ^ and u = R — p ^ 0, the 
Green's function Ge diverges as luj^**, or equivalently as 

(^\rA — Pa\'^ + \^B — Pb\'^) I but the second integral in 
Eq. l|E5p is convergent in u ~ 0: 



dpGB (R+-,R- -;p,p) [^rcg(p) - V'rcg(-R)] 



du 



u 



u 



dtp. 



reg 



du 



O \u 



u=0 



(E12) 



U= I dpGB(R+|,R-|;p,p 



(E6) 



Performing the Gaussian integration over p |68j | one finds 



U 



(27r) 



3/2 



dt- 



sin(t) -ftcos(i)|3/2^ 



(E7) 



where 



W 



IR 



r 12 
2I 



R 



2t 



1 



2 tan(t) X 



R 



R 



2t 



2sin(<) 



(E8) 

We have here introduced the shorthand notation x{t) — 
[1/i+l/ tan(i)]/2 and the phase factor (j){t), which equals 
exp(i7rn/2) for i„ < t < t„+i, where to = and 
ti = 2.029, ^2 = 4.913, . . . are the consecutive solutions 



(the first order term vanishes due to spherical symmetry). 
In this term therefore no divergence arises and we may 
set r ~ 0. The angular integrations can be performed 
analytically and wc obtain: 

/ dp Ge (R, R; p, p) [ipi-cg{p) - ipi-cgiR)] 

1 /"DO n 1 

with 

h (i|sm(t)|)"" (2»fs 

which is symmetric under the exchange of p and R. We 
have here introduced the radial wave function ^(-R) ~ 
Rtp{R) and the function y{t) = [l/t+ l/sin(i)]/2. 
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Writing 

i, (rA, vb) '■=° Vrcg(i?) (l - 7) + 0(1)' (E15) 

wc can cancel the divergent contribution on both sides of 
Eq. (jESp . and the remaining terms constitute the impUcit 
integral equation 



setting T = it/h in Eqs. HE2IE3IE4P and E = Zhwj^ in 
Eq. ljE4p : the derivation of O proceeds in an analogous 
way, but the resulting integral equation is much easier to 
solve since the integrand in both F\ and Fi become real 
damped functions with no finite time singularities. For 
an arbitrary value of the mass ratio a = mA/'niB^ they 
are given by: 



+ gj dpF2{R,p) [V'rcg(p) - |^\-cg(i?)] (E16) 

that needs to be solved numerically in order to determine 
i/'iog- The latter equation can be written in a symbolic, 
more compact form as 



V'r. 



I-gO 



aV'o 



(E17) 



where O is a symmetric integral operator, which is real 
in the limit k -^ 0. 

If one only aims at calculating strictly zero energy 
properties (i.e. Oeff, and not fk or rg), the treatment 
presented in this Appendix can be drastically simplified 



F^{R)=- 



(27r)3/2 7o 
exp(3T/2) 



dri -- 



,3/2 



+ 



[sh(T) + ^ch(r)] 



3/2 



exp —R 



, ch(T) + ^sh(T)-l - 

sh(r) + Jch(T) , 



and 



/•OO 

F2{R,p) =-a^/^ / dr 



(E18) 



exp[-(7?2+p2-)^] sh(2i?/9y) 
[rsh(T)exp(-T)]3/2 {2TTfy ' 

(E19) 



with2:(T) = a/2T+l/2th(r) 'AnAyir) = a/2T+l/2sh(r). 
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where fj, = m^ + m^ , fci, ^2 are the wavevectors 
of the two A atoms and fcs is the wavevector of the B 
atom. Whereas low enough values of fci , ^2 < kp may be 
chosen by adjusting the density of A atoms, the value of 
fca is fixed by the zero-energy oscillation of the trapped 
B particle, ft^fc|/2m_B ~ Tku. The fci = ^2 = limit 
of Eq. lE20t therefore results in a condition that looks 
difficult to satisfy in general for a^s ~ d: 



hcu- 



2mA 



h^ 



TUB -f 2mA 



2M«off 
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In the limit mA <^ rnB, however, as it is the case for 

the Li-Rb mixture, the vibrational energy in Eq. I|E21|I 

is multiplied by a small number ~ 2mA /m,B whereas 

the binding energy e increases as I/ttia: for aeS = d = 

Tv/kL and Vif = 50E^ , we thus find that the condition 

is satisfied for A = "^Li, B = ®''Rb, so that Eqs. (11,12) 

of [23 may be used. 

The T matrix of the matter wave is related to M by 

T(E + iO+) = (2^fiV"^A)Ei.J^"'Ldr.-)(rd. 

J. M. Ziman, Models of Disorder (Cambridge University 

Press, Cambridge, England, 1979), chapter 8. 

This results from the fact that the Hamiltonian TC is real 

in position space. Let \^) = G{E + irj)\ro). By definition 
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ofe, 

{E + ir^ ~H)iP{r)= 5 (r-ro), 

where 7i is a real differential operator. Taking the limit 
ry ^ O""" of the imaginary part of the latter equation leads 
to H0 — E(j>. Alternatively, using the spectral expansion 
of the Hamiltonian and of the resolvent with real eigen- 
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Fin as before. 

Requiring that the modulus square of the wavefunction 
decreases by four orders of magnitude over the radius of 
the medium (which is a minimal requirement to pretend 
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where (t>out is obtained from the far field extrapolation. So 
showing that this (jf>out is small amounts to showing that 
all the amplitudes ci^rn in the various output channels 
I, m are small. 

We also note that the intuitive fact that |77io(fc)|^ is min- 
imal in fc = fco leads to a real value for j3. 
Note that Hem'o has to be positive for consistency of 
the present 1ow-|7T1o| theory: with the reasoning lead- 
ing to Eq. IDlll . we get the approximation for the re- 
solvent {r\g{E + iO+)\ro) ~ -AF{r)F{ro)/mo{k). The 
linearization procedure in k around fco leads to mo{k) ~ 
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